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Abstract. We consider a model of active Brownian particles with velocity- alignment 
in two spatial dimensions with passive and active fluctuations. Hereby, active 
fluctuations refers to purely non-equilibrium stochastic forces correlated with the 
heading of an individual active particle. In the simplest case studied here, they 
are assumed as independent stochastic forces parallel (speed noise) and perpendicular 
(angular noise) to the velocity of the particle. On the other hand, passive fluctuations 
are defined by a noise vector independent of the direction of motion of a particle, and 
may account for example for thermal fluctuations. 

We derive a macroscopic description of the active Brownian particle gas with 
velocity-alignment interaction. Hereby, we start from the individual based description 
in terms of stochastic differential equations (Langcvin equations) and derive equations 
of motion for the coarse grained kinetic variables (density, velocity and temperature) 
via a moment expansion of the corresponding probability density function. 

We focus here in particular on the different impact of active and passive fluctuations 
on the onset of collective motion and show how active fluctuations in the active 
Brownian dynamics can change the phase-transition behaviour of the system. In 
particular, we show that active angular fluctuation lead to an earlier breakdown of 
collective motion and to emergence of a new bistable regime in the mean-field case. 



PACS numbers: 05.40.Ca, 05.40Jc, 87.10.Mn, 05.70.Ln, 05.20.Dd 
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1. Introduction 

The fascinating self-organisation phenomenon of collective motion in biology, such as 
exhibited by flocks of birds or schools of fish, has attracted scientists from very different 
backgrounds over the past decades. Whereas biologists focus on the evolutionary 
advantages and the biological function of moving in groups [IH3], physicist address 
rather the question about universal laws, governing the collective dynamics either from 
the theoretical point of view [H4T0] or, as done recently, by analysing certain model 
systems experimentally [TTHT3]. In general, physicist are inclined towards the analysis 
of simple, but generic, models of collective motion such as the minimal model proposed 
by Vicsek et al [4] . The Vicsek model can be considered as a non-equilibrium extension 
of the XY-model, where self-propelled particles move with a finite speed in continuous 
space and at discrete times align their velocities with the average velocity of other 
particles sensed in their local neighbourhood in the presence of noise. The Vicsek model 
exhibits a non-equilibrium phase transition leading to the onset of long range order at 
sufficiently low noise and high particle density. Detailed studies of the Vicsek model 
have led to an intense debate on the nature of the phase transition (continuous vs. 
discontinuous) [7H3C5UI6] . 

One of the main characteristics of the Vicsek model, as well as of many other 
individual-based models of collective motion, is the assumption of a constant speed of 
individuals. In general, this does not hold for the motion of real biological agents [17] . 
Furthermore, recent experimental realizations of swarming in non-living systems are 
based on rectification of fluctuation due to asymmetric friction and therefore imply a 
non- negligible role of speed fluctuations [TT|[T2"]. 

A possibility to introduce speed fluctuations into the dynamics of active particles 
is the addition of random increments from a well-defined probability distribution with 
zero mean to the speed of an individual in numerical simulations [T8UT9] . This approach 
is useful to analyze the robustness of the collective dynamics with respect to noise. 
However, in this kind of models the speed dynamics is always independent of the 
interactions between individuals. 

A different approach is based on the formulation of stochastic equations of motion 
for the positions r, and velocities Vi of N interacting individuals (i — 1, . . . , N): 



Here, m, is the mass of the focal individual, F{ is an effective force vector governing the 
deterministic dynamics, whereas r) i models an effective stochastic force, which account 
for the randomness in the motion of the individual. In general, both terms can depend 
on the position and velocities of all other individuals. In addition, the stochastic force 
/7j, must, by definition, explicitly depend on time, which is not necessarily the case for 
the deterministic force J^. 



dt 
dvi 



Fi + Vi 



(1.1a) 
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From a physical point of view, the velocity equation can be considered as 

a Langevin equation with actual physical forces on the right hand side. However, the 
full description of all physical forces determining the motion of an individual will be 
extremely complex and not very practical in most cases. Thus from a more general 
point of view, we can consider Ft and r\ i as effective forces responsible for accelerations 
of the focal individual in response to external stimuli and/or due to internal decision 
processes. The functional form of F^ and r) i can be set using reasonable modeling 
assumptions based on qualitative empirical observations [HElQ], or, if available, directly 
on experimental measurements of interactions between individuals [201 - 422] . 

This class of variable speed models not only overcomes the problems mentioned 
above, but also allows a straight forward mathematical analysis of the dynamical 
behaviour due to the usage of stochastic differential equations. We refer to active 
particles with self-propulsion, where both, the magnitude and the direction of the 
velocity are (stochastic) degrees of freedom, as active Brownian particles (ABP). This 
allows a clear distinction from so-called self-propelled particles (SPP), which usually 
refers to active particles moving with constant speed [U[23]. Originally, ABP referred 
to Brownian agents at far-from-equilibrium conditions with an internal energy depot, 
which can take-up energy from the environment and convert it into kinetic energy of 
motion [2H42T] . However, as non-equilibrium conditions and local energy conversion 
are characteristic features for all active particle systems, we suggest the above general 
definition of ABP. 

Recently, two of the authors considered non-interacting ABP with so-called passive 
and active fluctuations [28]. Passive fluctuations are assumed to be independent of the 
direction of motion, and may be motivated for example by a fluctuating environment. 
This could be e.g. thermal fluctuations, but in general also other non-equilibrium 
processes. Active fluctuations, on the other hand, are assumed to be correlated with 
the current direction of motion of the active particle. In the simplest case, they consist 
of independent stochastic processes parallel and perpendicular to the instantaneous 
direction of motion. Such fluctuations may be associated with fluctuations of the actual 
propulsion mechanism and are, per definition, a pure far-from-equilibrium phenomenon. 

The main focus of this work is the analysis of the impact of different fluctuation 
types on collective dynamics of interacting ABP. Here, we consider a generic, effective 
alignment interaction between individual particles, which represents a simple extension 
of the Vicsek model to the variable speed case. In particular, we will derive and 
analyze kinetic equations of ABP with velocity alignment and both passive and active 
fluctuations. Please note that already in the context of the Vicsek model, it was argued 
that different noise implementations may change the nature of the phase transition 
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2. Active Brownian particles with velocity-alignment 

2.1. Microscopic description 

We describe the active Brownian particle gas by a system of coupled stochastic 
differential equations. The dynamics of the i^ 1 point-like particle (i = 1,...,N) with 
mass 771, is determined by the following set of Langevin equations interpreted in the 
sense of Stratonovich: 

— =v i} (2.1a) 
■ 

m i~^f = -l( v i) v i + A* ( u Ai ~ vi) + »7 4 (t) . (2.16) 

The self-propulsion of active particles is described by a velocity dependent friction 
coefficient 7(1^), which can take both negative and positive values. Here, we use a 
piecewise linear function, the so-called Schienbein-Gruler friction [251130]: 

- 7(«i)t>j = -a ( 1 - t-^t ) Vi = a(v - \vi\)e iv . (2.2) 



\ V i\. 

The direction of motion of the 1^ particle is indicated by a unit vector e^ v (f) = Vi/\vi\. 
Particles moving with small speeds \v{\ < vq accelerate into their direction of motion 
whereas their velocity is damped for high speeds \vi\ > v leading to an average velocity 
unequal to zero. The relaxation time of the velocity dynamics is given by r Q = ma" 1 
in the deterministic case. 

The above friction function is apolar in the sense that the particles do not 
distinguish between their front and back. A corresponding polar friction function with 
linear velocity relaxation was discussed by Romanczuk and Schimansky-Geier [28] and 
was used previously in various models of active particles (see e.g. [31]). In the majority of 
cases where backward motion is negligible, e.g. in the case of moderate noise intensities, 
both models are equivalent. 

Due to the linear relaxation to the fixed point vq, friction function ( 12.2ft can be 
thought of as Taylor series of an arbitrary friction function around its stationary speed 
where 7(1*0) = 0. 

The interaction of particles is considered to be a linear velocity- alignment 
mechanism with characteristic time = mpT 1 . In (12.1 ty . the local mean velocity 
is denoted by 

UA i = jj~J2 v 'j ■ ( 2 - 3 ) 

Al 3=1 

The focal particle i aligns its velocity with the average velocity of all particles iV^ 
within a metric interaction region Ai around its current position. In general, can 
have different geometries, however, a common choice in two spatial dimensions is a disc- 
like interaction zone with radius e. The velocity- alignment strength is denoted by jl. 
The summation in ( 12. 3 p includes the i^ particle itself so that in case of no neighbours 
(Na, = 1) the alignment force vanishes, consistently. If the local mean velocity ( 12. 3p 



Active Brownian particles with velocity- alignment and active fluctuations 



5 



equals zero, the alignment interaction reduces to Stokes friction i)i oc — JxVi decelerating 
the particles ("social friction"). 

In the limiting case of constant speeds \vi\ ~ vq, the alignment interaction is 
related to the Kuramoto model [32] describing synchronization phenomena of coupled 
oscillators. In this case, the polar angle defining the direction of motion corresponds 
to the phase of the oscillator [53^, 134]. Peruani et al derived a mean-field theory for self 
propelled particles with constants speed and Kuramoto interaction [35J. 

The above-mentioned velocity-alignment interaction has already been used to 
described collective dynamics, e.g. by Czirok et al [36] in order to describe the formation 
of complex bacterial colonies as well as by Niwa [37] to characterize the dynamics of fish 
schools based on Langevin equations. 

In ( 12.161 ). rj^t) denotes a stochastic force, introduced in [28], contributing three 
terms to the velocity dynamics. The vector points perpendicular to the direction of 
motion and has unit length one. 



The first expression on the right hand side denotes so called passive or external noise 
acting in the same way on both components of the velocity vector with the intensity 
d>E- It might be associated for example with ordinary Brownian fluctuations. Thus, 
in principle it can model essentially a equilibrium phenomenon linked to thermal 
fluctuations. In the present context, however, it may also be of non-equilibrium origin. 
We refer to this type of noise as passive, as it does not "carry" any informations on the 
propulsion of a particle. 

We refer to the remaining two terms together as active or internal noise. It is 
motivated by the assumption, that this internal fluctuations are directly linked to 
fluctuations of the driving "motor" or internal decision processes of the Brownian agent, 
i.e. animals. These active fluctuations are subdivided in two parts: speed noise, acting 
along the direction of motion with intensity d v and pure angular fluctuations, changing 
the direction of motion of the particle with intensity d v without affecting the particles 
speed. In the case of constant speed it is not possible to distinguish active and passive 
fluctuations as the angle defining the direction of motion is the only degree of freedom. 
Figure [T] illustrates the different effects of active and passive fluctuations on the particles 
velocities. 

The random processes £i tV (t) are assumed to by Gaussian white noise with zero mean 
and vanishing correlations: = 0, (€i, v (t)€j,p(f)) = 5ij5 Ufl 5(t - t') e {1, N}, 

v, /x e {x, y, v, <f}). Due to the non-equilibrium nature of the self-propelled dynamics and 
active fluctuations, the fluctuation-dissipation theorem does not hold here in general. 
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Figure 1. Visualization of passive (a) and active fluctuations (b). The dashed line 
shows the trajectory of the i th particles. Speed noise changes the speed of the particle 
only whereas angular noise acts on the direction of motion. On the other hand, the 
components of the velocity vector are changed by passive fluctuations in the same way. 



2.2. Natural units and diraensionless parameters 

First of all, we identify the relevant dimensionless parameters and characteristic length 
and time scales in ( 12. II) : 

(2.5) 



L = 


~mvo~ 


T = 


'171' 




. a . 




. a . 



Apparently, the velocity scale is defined by Vq. In this section, physical values in SI units 
are indicated by a prime. In order to reduce the number of parameters, the following 
transformation to dimensionless quantities r, v and t is used: 



mv 



t' 



in 



t 



v 



v v 



a a 

The Gaussian processes £(t') are transformed in the same manner: 



m 



(2.6) 



(2.7) 



Using these transformations, one can rewrite the microscopic Langevin dynamics (12.11) 
of the i** 1 particle in dimensionless quantities as follows: 

dt 



v, 



(2.8 a) 



= e ijV -Vi + fi [u Ai - Vi) + \/lD E + v 7 2 A, ii,v ei, v + y/2D v £ i>v e ijlp . (2.85) 



In (12.81) . the dimensionless velocity- alignment strength \x and noise intensities D v were 
introduced [y G {E, v,ip}): 

a 



D„ 



mav, 



(2.9) 



o 



Hence, the dynamical behaviour is governed by these four dimensionless parameters. 
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3. Hydrodynamic theory of active Brownian particles 



A hydrodynamic description based on symmetry properties of flocks was first proposed 
by Toner and Tu (TT) [HE]. They proposed a general continuum model for collective 
motion without considering a particular microscopic dynamics and determined the 
scaling exponents for density and velocity fluctuations. Hydrodynamic equations for 
SPP with aligning collisions were derived by Bertin et al in [38j[39] using a Boltzmann 
approach, which allows to link microscopic model parameters to the parameters of a 
coarse-grained theory of TT-type. Bertin et al assumed a small collective velocity, 
which implies the validity of the theory only close to the transition point. Recently, Ihle 
derived a kinematic description of the Vicsek model accounting also for multiparticle 
interactions [ID]. Related hydrodynamic descriptions of self-propelled rods were also 
studied by Marchetti and coworkers pTTlEO] . 

In order to predict the macroscopic behaviour of the active Brownian particle gas, 
a hydrodynamic description is derived directly from the microscopic Langevin dynamics 
in the following. We use a moment expansion of the corresponding probability density 
function (see also Appendix A for further details). The approach used in this work 
goes back to Riethmiiller et al , who studied gravity-driven granular flow [131 . It was 



successfully used to analyze the mean-field dynamics of ABP with velocity alignment in 
one dimension [H] and in two dimensions with passive fluctuations only |45j . 



Let p (rj, Vj, t) d 2 Vjd 2 Vj be the probability to find the j^ 1 particle at time t in the 



velocity interval [vj, Vj + dvj] and within the spatial region [rj, rj + drj]. The dynamics 
of the one-particle probability density function p (rj, Vj, t) can easily be derived from the 
dynamics of the joint probability density function P = P (r±, r 2 , r^, v±, v 2 , Vn, t) 
by integrating over all particle positions and velocities except the j . Due to the fact 
that any j could be chosen, the particle index is omitted henceforward. The dynamics 
of the single-particle probability density function p(r,v,t) reads 
dp 



dt 



+ 



+ 



V r (vp) - V v 
2 



v 

\v\ 



V 



dv 2 
d 2 



D v v 2 x + D v vl 



P 



dv x dv,. 



\v\ 

'VxVy 
I |2 

\v\ 



p + H (u A -v)p-D H 
d 2 fD v v 2 y + D v v 2 x 



v 



\v\ 



V 



dv 2 



P 



17 



{D v -D v )p 



&v(D E p) . 



Moments of the velocity components are defined by 



(vTv n 



J x y 



) 



P(r,t) 



d 2 v v^p(r,v,t). 



The spatial density of particles corresponds directly to 



P{r,t) 



J d 2 v p(r, v, t) . 



Multiplying fl 3 . 2 [) with the density and taking the temporal derivative yields 



(3.1) 



(3.2) 



(3.3) 



(3.4) 
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We may rewrite the integral on the right-hand side of (13.41) by inserting the Fokker- 
Planck equation (13.11) and using the definition (13.21) to obtain a general equation of 
motion for the M^ m?l ^-moment of the probability density function: 

— (o(v m v n )) = - — (p(v m+1 v n )) - — ( 
dt ^ x y/) dx {PK x y/) dy [ 



+ (m + n)p 



1 i ?^ 

u x u y 

\v\ 



1 ) 7 )^ 

u x u y 



\V\ 



+ fip [m u A , x (v™ l Vy) + n u A)V (v™v% x ) - (m + n) 



{«)] 



+ m(m — l)p 
+ n(n — l)p 



? .m-2 ? .n+2 
u x u y 



\V\ 

m+2 n-2 
u x u y 

\v\ 



D v 



7 )^^7 t^ 1 

u x u y 



\V\ 

u x u y 



\V\ 



f D E «~ 2 <> 
D E <«- 2 > 



+ 2mnp ( 



1 J 



\ \v\ 



(3.5) 



Up to now no explicit approximations were made except the mean field assumption 
(N Ai » 1). 

In order to derive hydrodynamic equations from (13.51) . we rewrite the individual 
velocity v as a sum of a mean-field velocity u(r,t) at the position of the focal particle 
plus some deviation dv(r,t): 

v — u(r,t) + Sv(r,t) . (3.6) 

We assume that the expectation value of 6v% (r,t) vanishes for odd k: 

(5v k Xyy ) = for k = 1,3,5,... . (3.7) 

This corresponds to the assumption of a symmetric p(r,v,t) with respect to u(r,t) in 
the velocity space. Furthermore, we assume 



({5v x ) 2 Sv y ) = ({5vy) 2 5v x ) = 



(3. 



In addition to the mean-field velocity, we introduce the covariance matrix as a measure 
of velocity fluctuations around the mean-field velocity: 

{(Sv x ) 2 ) (5v x 5v y ) 



(5v y 8v x ) ((8vy) ) 



)- 


f T 

- 1 XX 


T X y \ 






Tyy I 



(3.9) 



In the following we will refer to T as the "temperature tensor" . 

In the lowest order (m = n = 0), one obtains from (I3.5P the continuity equation for 
the local density 
dp 
di 

reflecting the conservation of the particle number N. 

To obtain evolution equations for the mean-field velocity and the temperature 
tensor, we need to evaluate the expectation values of fractional expressions of the form 



V r (pit) = 



(3.10) 
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Jy/\v\ k in (13. 5|) . Here, we use the approximation 



1 ?^^7 

u x u y 



\V\ 



yk 



(3-H] 



with V = V(u,T) = \J\u\ 2 + tr(T) . This approximation provides naturally a closure 
of the hierarchy of moment equations with m + n < 2. 

From equation ( 13.41) for the first order moments with m — 1, n — and m = 0, 
n = 1 combined with (13. lOj) . we can derive the evolution equation for the mean- field 
velocity field to 



du 

— + liV, u 

at 



1 

1 

V 



u 



v 2 



u 



fi (u A -u)-- { 

p 



V T ( P T) 



The above equation (13 . 1 2[) is a partial integro-differential equation since 



u A (r,t) = A" 1 [ dV u{r',t). 

J A 



(3.12) 



(3.13) 



Csahok and Czirok [16] have suggested a hydrodynamic velocity equation for collective 
motion of bacteria with similar active propulsion and interaction terms as derived in 
Eq. (13.121) . We should note that in general, for a finite numbers of neighbours u e , will 
be a fluctuating quantity which is neglected here by assuming N t 0. 

Finally, we can combine ( 13 . 1 0[) and ( 13. 1 2[) together with the three second-order 
equations obtained from ( 13.51) for m + n = 2 ((m, n) G {(2, 0), (0, 2), (1, 1)}). This yields 
the equation for the temperature tensor 



dT 

~dt 



{uV r )f 



1 

V 



1 

+ 2 



\T-fiT + D E I 
(TV)u T + | 



A 



V 2 9 



with / being the identity matrix and 



T 



[uu 



^x^y ^xy 



U X U 
U 



2 +T 



T 

J-xy 



!)!) 



u 



y 1 



T 
± yy 



T 

- 1 XX 



2T, 



2T 



< + T x 



T 



(3.14) 



(3.15a) 



(3.156) 



yy 



From ( I3.12p and (I3.14p one can already identify the general impact of the different 
fluctuation types: The active angular fluctuations D v introduce an additional effective 
friction which affects both the velocity and temperature dynamics, whereas the active 
velocity fluctuations D v and passive fluctuations D E drive the temperature equation. 

The kinetic theory contains the important limiting case of ordinary Brownian 
motion, which corresponds to /i, D v , D v v equal to zero. In this special case, we recover 
the fluctuation-dissipation relation known from equilibrium statistical mechanics. 
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(4.1) 



4. Spatially homogeneous, steady-state solutions 

Under the assumption of a spatially homogeneous system, we may drop all spatial 
derivatives in (13.101) . (I3.12p and (13.141) . Furthermore, in the steady-state case, we 
may choose the coordinate system such that u = (|u|,0) is parallel to the x-axis 
of the laboratory frame of reference. In this case, we assume that the average 
deviations parallel and perpendicular to the direction of motion become independent 
(Sv x 5v y ) = (8v\\6v±) = (6v\\) (Sv±), i.e. the temperature tensor assumes a diagonal 
form 

<(^||) 2 > <*«_!_> \ / T, 

(5v x ) (5v\\) ((5v ± ) 2 ) J [ T ± 

where Ty is the mean squared velocity deviation parallel to u, whereas Tj_ is the mean 
squared deviation perpendicular to u. We define T = tr (T) as the (scalar) kinetic 
temperature. The simplifying assumption of independent average fluctuations parallel 
and perpendicular to the average velocity is justified by the broken symmetry of the 
collective translational steady-state [6]. Static correlations in the macroscopic mean- 
field velocity fluctuations would imply collective rotational motion (see e.g. [28]), which 
are not stable in our model. 

In the spatially homogeneous and stationary case, we obtain the following set of 
homogeneous equations. 
d|u| 
~df 



0= l(\u\ 2 + T) > -l-D^duf + T) )\u\ (4.2a) 



0= ({\u\* + T)-t-l-n)T l+ D< ' 



2 dt \" 1 rj n ■ ^ y\ u f + T 

Dv[ w^) +De {A - 2b) 



+ D v (—^—)+D E . (4.2c) 
\|u| +TJ 

The corresponding stationary solutions of the ordinary differential equations (14.2 a§ - 
(14.2c| can in principle be found but are quite complex and hard to analyze. A simpler 
equation for the temperature can be found by adding up (14.2 £>|) and (14.2 cp . which yields 



d\u 



dt 

ldT 



0= ((|u| 2 + T) 5 - 1 - D v (|u| 2 + T) l ^\u\ , (4.3a) 

^ = 0= {[\\u\ 2 + Ty L2 -l-^T + D Lp \u\ 2 {\u\ 2 + Ty 1 + D v + 2D E . (4.36) 

In the following, we discuss the time-independent stationary solutions of (14.31) which 
correspond to fixed points of the dynamical system. One pair of fixed points reads 

Mi = 0, (4.4a) 
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T} ±} = 2D °^ V + ri~ 2 (l ± Vl + 4VD E + D v )(l + n) ) . (4.46) 

However, only one of these fixed points is physically relevant. Solution (^\u\ 1 , t/ will 
be neglected because it predicts unphysical behaviour for vanishing noise intensities. 
In this case, both \u\ 1 and t/ - -* are equal to zero. This corresponds to all particles 
being at rest, which is not possible for permanently self-driven particles considered 
here. Solution y\u\ 1 , T"/ + M with vanishing mean speed and finite fluctuations is called 
disordered state. It is stable above a critical noise intensity and below a critical coupling 
strength, respectively. The temperature in the disordered state is independent of the 
intensity of angular fluctuations D v . 

Another two solutions with nonzero mean speed are given by the following 
expressions: 



l»f = ^(l±V^45;)-C y - gE±^±jgg , (4.5a) 
T 2 = ^ + *>. + M« (4. 56 ) 

We refer to these fixed points as ordered state. For vanishing noise intensities 
(D E = D v = Dp = 0), one solution reads — 1 an d T 2 = 0. It corresponds to 
all particles moving perfectly aligned in the same direction with stationary speed of 
individual particles \vi\ — 1. In general, the solution u {+) is stable for sufficiently weak 
noise, or strong alignment, respectively. The second ordered solution u 2 ^ is always 
unstable. 

In ( 14. 31) . the intensities of passive noise D E and active speed noise D v enter the 
equation with the same sign. The theory predicts a qualitatively identical influence of 
passive fluctuations and speed noise on the mean-field behaviour of the system. 

For weak interaction \l < 1, the theory always predicts a continuous transition of the 
mean velocity from the ordered to the disordered state (figure [2^,b). In the case of strong 
interaction \x > 1 and not too large passive fluctuations, a new regime emerges. In this 
regime, both the ordered and the disordered state are stable macroscopic configurations 
(figure [2b, d). The theory predicts that, depending on the initial condition, the system 
relaxes into one of these stable macroscopic states. Therefore, a discontinuous transition 
from collective motion to the disordered state for increasing angular fluctuations is 
expected. In figure El we show the phase diagram with respect to external noise D E 
and active angular fluctuations D v for /x = 10 and D v = 0. In the bistable regime, if 
the mean velocity ua = u is low, the velocity-alignment force reduces to Stokes friction 
(social friction) which dominates the individual particle dynamics and as a consequence 
leads to a stable disordered state. Surprisingly, the transition point D VtC = 0.25 does 
not depend on the coupling strength fj,. From the Fokker-Planck equation ( 13.11) one can 
see, that D v > leads to an effective drift against the collective direction of motion of 
the particles. Subsequently, the individual particle velocities decrease and above D V}C 
no self-driven motion along the consensus direction motion is possible. 
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Figure 2. Comparison of numerical Langevin simulations of (|2.8|) (points) and 
theoretical prediction (|4.3[) (lines) in the case of vanishing speed noise and passive noise 
(De = D v = 0). The dashed lines indicates unstable states whereas the solid lines 
represent stable configurations. Results of simulations initially started in the disordered 
state are indicated by (black) squares. Stars (red) represent results of simulations which 
were performed using the ordered state as initial condition. The simulations were 
performed considering global interaction Aj = L x L, periodic boundary conditions 
and N = 10 4 particles. Parameter: \i = 0.1 (a)-(b), fi = 10 (c)-(d), numerical time 
step At = 1 • 10~ 4 . 

Please note that the mean speed in the limiting case of strong interaction fi — > oo 
still depends on the intensity of the angular fluctuations. 

lim \u\^ = \ \ (l ± y/l-4D v ) - D v (4.6) 

A»— too V 2 \ / 

The analytical results (I4.4p and (14. 5 p are compared with numerical simulations of 
the Langevin dynamics (12.81) . We used the stochastic Heun-Scheme for the numerical 
intergration, which converges in quadratic mean to the solution of Langevin equations 
interpreted in the sense of Stratonovich ^7]. The simulations were performed with 
N = 10 4 particles considering a square spatial area of length L and periodic boundary 
conditions. To ensure a spatially homogeneous system, global interaction was used 
(Ai = L x L). In order to detect the bistable region, both the ordered and disordered 
state were used as initial condition. Results of the numerical simulations in comparison 
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Figure 3. Phase diagram in the case of global interaction for D v = 0. Three different 
regimes can be distinguished: only ordered solution stable (A), only disordered solution 
stable (B), and a bistable regime (C). The bistable regime, where both the ordered 
and disordered states are stable, exists only for finite active angular noise intensities 
and strong interaction u > 1, e.g. for De = 0: fj,(u + 1)~ 2 < D v < 0.25. It does 
not exist for the case of passive fluctuations only [45]. The solid lines represent the 
borders of the different regimes as predicted by the theory. The symbols represent the 
results of numerical simulations with TV = 10 4 particles. At each point in the phase 
space two simulation runs were performed with different initial conditions (ordered 
and disordered) until the system reached a stationary state. The circles represent the 
ordered stationary result, whereas the crosses show the results with a disordered final 
state. The dashed line shows the critical active angular noise D VyC — 0.25. Numerical 
time step in all simulations At = 5 • 10 -3 . 



with analytical predictions of the mean-field theory are shown in figures [2] and [3j 

In the absence of active fluctuations, the theory is mainly confirmed by the 
numerical simulations. This case was discussed in detail in jl5]. If the intensity of 
the (active) speed noise is small in comparison to intensities of the active angular 
and passive fluctuations, the theoretical predictions are in good agreement with the 
numerical results. 

Deviations from the theory are due to the approximations used, especially (13. lip . 
This approximation implies that the dynamics of the local velocity field follows closely 
the individual dynamics or more specifically: the decrease of the mean speed is related 
to a decrease of the individual particles velocities. Thus, the agreement of theory and 
simulations is better for a strong velocity- alignment, whereas small deviations are found 
in the case of fixed particles speeds (// <C 1). 

According to the phase diagram in figure [3], it is possible to move from the bistable 
regime (B) into the only ordered regime by increasing the external noise strength De- 
This implies the rather counter-intuitive effect of "generating" order from disorder by 
increasing noise. If we consider the system to be initially in the disordered state in the 
bistable regime, an increase in De will eventually lead to a destabilization of the current 



Active Brownian particles with velocity- alignment and active fluctuations 



14 



s" 



0.75 

0.5 

0.25 


1 

0.75 
0.5 
0.25 




1 1 1 1 


rz. ~ 


P 
o 

0.2 K- 

CD 






CB 
err 






:eng 




mean speed u 


err 




order parameter $ 


0.02 J 3 




— noise De 



mean speed u 

order parameter $ 

— noise De 



50 



100 150 
time t 



200 




- 0.2 



0.02 



250 



o 



o 

crq 

b 

5] 



Figure 4. Switching of the globally coupled system from the disordered state to the 
ordered state with increasing noise intensity for step like noise change (top) and a linear 
switching protocol (bottom). The system is initialized in the disordered state with 
a noise intensity of the passive fluctuations De = 0.02 corresponding to the bistable 
regime. At time t = 100, the noise intensity is increased either due to an instantaneous 
step-like switch or in finite time via a linear protocol to the final value De = 0.2 which 
corresponds to the ordered regime. In addition to the mean speed |tt| = u, we show 



also the directional order parameter $ = 
parameters: D v = 0.2, D v = 0, N = 10 4 , /z 



E- =1 

= 10. 



/N with ei V = Vij 



Other 



state and to a transition from disordered to ordered solution. The actual transition 
depends on the protocol D^if) as shown in figure HI However, if the noise increase is 
not to strong, the system will eventually relax to the ordered state of collective motion. 

We point out that the found bistable behaviour differs from the previously reported 
bistabilities in collective motion (see e.g. [441148] ) . which were either a consequence of a 
non-linear friction function in one-dimensional dynamics jB] or were caused by Morse- 
type interactions of active particles in three dimensions [48]. Noise in the cited cases 
did not create the coexistence of the two phases but merely induced switching between 
the two states with switching times which appeared to be exponentially distributed. 

5. Discussion &z, Outlook 

In this work, we focused on the impact of the different fluctuation types on the 
collective dynamics of active Brownian particles with velocity alignment. We derived a 
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hydrodynamic theory for the coarse-grained observables spatial density p, velocity-field 
u and the temperature tensor T in the presence of both fluctuation types. However, 
approximations are necessary in order to obtain a closed set of equations, which may lead 
to deviations of theoretical predictions from the numerical simulations. For simplicity, 
we restricted ourselves to stationary, spatially homogeneous systems where analytical 
results can be obtained. This case corresponds either to global coupling (interaction 
length comparable to the system size) or to local dynamics of the system at length- 
scales where gradients in the macroscopic variables can be neglected. Our analysis 
shows that the introduction of active fluctuations has a striking influence on the mean- 
field behaviour of active particles with velocity-alignment and variable speed: In the 
case of strong interaction \i > 1, a new bistable regime emerges that cannot be found 
for passive noise only. The corresponding transition from the ordered to the disordered 
state is discontinuous for increasing angular noise, whereas it is always continuous in 
the mean- field for purely passive noise |45j . Our analytical results agree well with 
the simulation results for strong alignment ji ^> 1 and low active speed fluctuations 
D v which is related to the approximation used for the expectation values of rational 
expressions appearing in the general moment equation and the corresponding implicit 
moment closure. 

The onset of order in the spatially homogeneous system is related to the 
synchronization of coupled oscillators j33j[3l]. If we identify the modulus of the velocity 
with the amplitude and the direction of motion with the phase of an oscillator, the 
active particle system with velocity- alignment corresponds to a system of coupled noisy 
oscillators with the same eigenfreqency in a co-rotating frame of reference. In the limit 
of small coupling /i C 1, the critical noise strengths for the transition from order to 
disorder have to be small as well (D E , D v , D v <C 1). Thus, the fluctuation in the modulus 
of the velocity at the critical point will be negligible. The system reduces effectively to 
coupled identical phase-oscillators which is a special case of the well studied Kuramoto- 
model [32]. In this context, our results for strong alignment can be linked to amplitude 
death in system of strongly coupled limit-cycle oscillators [19]. The effective social 
friction in the disordered state due to the velocity alignment slows down the individual 
particles and stabilizes the disordered regime. 

We focused here on the case of global interactions, however the hydrodynamic 
theory provides the basis of the stability analysis of the spatial inhomogeneous solution. 
It is known from the Vicsek model that density fluctuations occur in the case of 
local velocity-alignment interaction [TJ. Recently, a number of studies analyzed the 
coarse-grained dynamics of Vicsek-type models showing the instability of the ordered 
homogeneous state of collective motion close to the transition point [38VH01I50] . Similar 
density inhomogeneities can be observed for active Brownian particles with local 
coupling and passive noise for sufficiently small interaction regions (Ai <C L x L) [45J. 
Although our results for the homogeneous system in the presence of active fluctuations 
cannot be directly generalized for the case of local interactions, they are nevertheless 
relevant for it, as they provide predictions on local dynamics in high density regions. In 
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the regime where only the ordered solution is stable (e.g. for purely passive fluctuations), 
regions of high density are correlated with regions of high order. Thus, they correspond 
to clusters or bands performing collective translational motion. This is also the case for 
the classical Vicsek model [7] . However, the situation changes drastically in the bistable 
regime. Here, high density structures may be both ordered (mobile) or disordered and 
stationary. 

In figure we show an example of a simulation for local coupling. We assume that 
a focal particle interacts with its neighbours within a metric range e = 1 (L = 100). 
The parameters are chosen to correspond to the mean-field bistable regime (/i = 10, 
Dtp = 0.2, D v = 0.0, De = 0.01). After a short initial transient, the system reaches 
a metastable state of multiple stationary clusters at t ~ 1000, with few small mobile 
clusters in between (figure [5A). At time t ~ 11000 a sufficiently large mobile cluster 
forms due to fluctuations. It is able to "break up" the stationary high density aggregates 
and we observe a sudden onset of large scale collective motion (figure [5]A.). At time 
t « 20000, we observe again a spontaneous formation of a stationary (disordered) cluster 
due to a collision of clusters moving locally in different directions at x ~ 70, y ~ 30 
(figure [5]C) . This cluster represents effectively an absorbing obstacle for the remaining 
mobile structures and leads eventually to a break-down of collective motion (figure [5p). 
However, at the end of the simulation run at t m 40000 we observe again the onset 
of collective motion (figure EE) • Clearly, the bistable local dynamics predicted by the 
mean-field theory, leads to complex spatial and temporal dynamics for local coupling. 
Please note, that the time scales governing the large scale dynamics are in the order 
~ 10 4 . This is five orders of magnitude larger then the interaction time scale fi^ 1 =0.1 
which governs the microscopic dynamics. This large spread in the relevant time scales 
makes this system difficult to analyze by standard numerical methods and may require 
novel computational techniques, as for example the equation- free approach [5lTf53] . 

In the context of local coupling, we should mention recent studies on collective 
dynamics in different variants of the Vicsek model [5UE5] and in a run-and-tumble 
model [56]. In these models the velocity of individuals was assumed to be directly a 
function of the local density of individuals [55JE6] , or of the local order parameter [54]. 
These models show qualitatively similar phase-separation dynamics, resembling the 
behaviour observed here for local coupling. However, in our model we do not explicitly 
include a dependency of the individual speed on the mean-field velocity and density, 
but it emerges naturally due to the interplay of the velocity-alignment force and the 
variable speed of individuals. 

In conclusion, the ability of active Brownian particles to adjust their speed in 
response to interactions with the environment is certainly a more realistic, in comparison 
to the constant speed restriction in self-propelled particle models. However, this 
extension is far from being trivial. The mean-field approximation gives important 
insights in the emergent collective behaviour of the system and demonstrates clearly 
how the type of fluctuations can have a dramatic impact on the large scale collective 
dynamics. 
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time t 

Figure 5. Active Brownian particles with local coupling (disc-like interaction zone 
Ai with interaction radius e = 1, L = 100, N = 10000) in the bistable regime. (A-E): 
Spatial snapshots at t = 10000 (A), t = 15500 (B), t = 20000 (C), t = 30000 (D) and 
t = 39900 (E). Each particle is shown by its velocity vector, with the color indicating 
the direction of motion according to (F). (G): Time evolution of the total mean velocity 
U(t) = \ J2 i = l N vA of the corresponding simulation run with arrows indicating the 
snapshot times. The initial condition was a random distribution of particles and their 
directions of motion (disordered state). Other parameters: D v = 0.2, D v = 0.0, 
De = 0.01, n = 10, numerical time step At = 5 • 10~ 3 . 



Appendix A. Derivation of the moment equation 



The general evolution equation (13.51) for the moments M^ mn ^ is derived in this section. 
Therefore it is exploited, that the dynamics of the active Brownian particle gas, which 
is determined by four stochastic differential equations (12.81) . can equally be described by 
the corresponding Fokker-Planck equation, which describes the dynamics of the total 
probability density function P (t*i, rjv, V\, ujy, t). The Stratonovich interpretation 
of the stochastic differential equation is used, which yields: 
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(A.l) 

(A.2) 



hence (1A.1I) is a nonlinear Fokker-Planck equation. The interaction region of the i^ 
particle is denoted by A±. 

In the main text, the one-particle probability density function p(r,v,t) is 
introduced. Its dynamics can easily be derived from (lA.ip by integrating over all particle 
positions and velocities except one, arbitrary chosen particle as well as using the mean 
field approximation 

j A d 2 r J d 2 v vp(r, v, t) 



u A (r,t) 



(A.3) 



J A d 2 r J d 2 v p(r,v,t) 

This assumption certainly holds for high particle densities and leads to an effective 
single-particle description: 
dp 
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Multiplying the moments of the velocity-components 
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' dt 

We may rewrite the integral on the right hand side of (1A.7|) by inserting the Fokker- 
Planck equation (IA.40 
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Integrals containing a derivative with respect to v x and v y , respectively, are integrated 
by parts and rewritten using the moment definition (IA.5[) . Furthermore it is assumed 
that p(r,v,t) decreases fast enough, so that 
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holds. In this way, partial differential equations for the moments (macroscopic variables) 
of the probability density function are derived directly from the microscopic Langevin 
description (12. ip . The general evolution equation for the moment M^ mn ^ reads 
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It is shown in section [3] how a closed set of hydrodynamic equations up to the second 
order can be derived from ( 1A.10|) . 



Appendix B. Stochastic Heun method 



A system of stochastic differential equations with multiplicative Gaussian white noise 

«&(*)> = o, (fc(f)&(f)> = M(*-0), 

x = f(x,t) + g(x,t)$(t), (B.l) 

can numerically be integrated applying the stochastic Heun method. This integration 
scheme converges in quadratic mean to the solution of Langevin equations interpreted 
in the sense of Stratonovich jl?] and reads as follows. 



x 



(t + At) = Xi (t) + fi(x(t),t)At + 9ij(x(t), t)aj{t)y/Ai 



(B.2a) 
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Xi (t + At) = Xi (t) + - [fi(x(t),t) + fi(x{t + At),t + At)] At 



^ m 



2 

i=i 

For the integration, Gaussian random variables with zero mean and variance one are 
used. 

= (B.3a) 

(<7i(*K-(*)) = 5ij (B.36) 
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